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Abstract 

We consider the single-inclusive minijet cross section in pA at forward rapidity 
within the Color Glass Condensate model of high energy collisions. We show 
that the nucleus appears black to the incident quarks except for very large 
impact parameters. A markedly flatter pt distribution as compared to QCD 
in the dilute perturbative limit is predicted for transverse momenta about the 
saturation scale, which could be as large as ~ 10 GeV^ for a gold nucleus 
boosted to rapidity ~ 10 (as at the BNL-RHIC). 

QCD correctly predicted logarithmic violations of Bjorken scaling in Deep- Inelastic 
electron-proton scattering at very large Q^, i.e. at very short distances Asymp- 
totic freedom provides the theoretical basis for the successfull applications of per- 
turbative QCD to hard scattering, short distance phenomena. However, the region 
of QCD phase space where the field strengths are strong is largely unexplored. This 
is where one expects that cross sections become comparable to geometric sizes of 
hadrons and nuclei (the "black limit") and where the unitarity limit is reached. A 
perturbative QCD based mechanism for unitarization of cross sections is provided 
by gluon saturation effects 0, 0]. A semi-classical approach to gluon saturation and 
QCD in the high energy limit (small x) was developed in ^ |^, 0| and applied to 
high energy heavy ion collisions at RHIC P, |^, [1^, |Tl]]. 

An ideal environment to study gluon saturation and unitarization effects is pro- 
vided by large nuclei since high gluon density effects are expected to be amplified by 
a factor of A^^^ due to the Lorentz contraction of the nucleus. The scale associated 
with the high gluon density, the saturation scale Qg, grows with energy and A and 
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decreases with increasing impact parameter. The Relativistic Heavy-Ion CoUider 
(RHIC) at BNL will soon allow experimental study of proton-gold or deuteron-gold 
collisions at a center-of-mass energy of ^/i ~ 200 — 300 GeV. We suggest that the 
saturation regime of QCD can be probed at RHIC by measuring the inclusive cross 



section in p + Au (or d + Au) collisions (in this regard, see also ||Tl| , In par- 

ticular, in the forward region, i.e. close to the rapidity of the proton beam, the 
saturation scale Qs can become quite large due to renormalization group evolution 
in rapidity 0. Thus, we predict significant modifications of the pt distribution of 
produced pions relative to leading twist perturbation theory at transverse momenta 
as large as several GeV. A modification of the longitudinal distribution of leading 
hadrons produced by electrons scattering inelastically from a black target has been 
predicted previously |13[. Here, we focus on the transverse distribution in the for- 
ward region from p + A scattering, which will be analyzed experimentally in the 
near future at RHIC. 

At large rapidity, we consider the quark-nucleus elastic and total scattering cross 
sections. (In turn, towards midrapidity gluon production becomes the dominant 
contribution to the cross section in the Color Glass Condensate model We 
argue that the total quark-nucleus scattering cross section may be related to the 
single inclusive hadron (jet) cross section in proton-nucleus collisions by using the 
colinear factorization theorem on the proton side. Let p^ {q^) be the momentum 
of the incoming (outgoing) quark. We assume the quark is moving along the left 
branch of the light cone such that p~ ^ p^ = p1/2p~ . The starting point is the 
scattering amplitude (for brevity, we do not write polarization indices explicitly) 

{(l{(l)out\q{p)in) = {out\bout{q)blr^{p)\in) , (1) 
which, using the LSZ formalism [|14| can be written as 

{out\KUq)blmn) = -^Jd''xd''ye-'^P^-'^yMQ)[^0y-m] 

{out\T'ilj{y)ijj{x)\in)[—i 0^ —m]u{p) (2) 

where m is the quark mass and Z2 is the fermion wave function renormalization 
factor. u{p) is the quark spinor with momentum p. The fermion propagator Gi;' in 
the background of the classical color field is 

{out\Tip{y)ijj{x)\in) = —i{out\in)Gp(y,x) . (3) 

The amplitude then becomes 



{cii(l)out\qip)in) = —{out\in) I d xd ye 
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(4) 



In momentum space, the fermion propagator Gp can be written as H, nB, 16 



GF{q,p) = (27r)V(g (p) - zgG%{q) 



(2vr) 



-Mk)GF{q + Kp) 



(5) 



where 4 = ^'^T/i is the classical background color field, and is the free fermion 
propagator. It is useful to define the interaction part of the fermion propagator from 



as 



Gp{q,p) = {2nyS\q-p)GUp) + G'AqMq,p)GUp) • 
Substituting (|) into the amplitude @) leads to 

{(l{<l)out\q{p)in) = u{q)T{q,p)u{p) . 



(6) 
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Here, we have set Z2 = 1 and {out\in) = 1 since we are working to leading order 
in as and our background field is time independent. This is a very simple relation 
between the amplitude for scattering of a quark or anti-quark from the Color Glass 
Condensate and the quark propagator in the background color field of the nucleus. 

The explicit form of the quark propagator in the background of a classical color 
field was calculated in |16[. The interaction part, as defined in (P) is given by 

PI 



r{q,p) = (27r)5(p - g )7 / dht 



Vizt] 



Jiqt-pt)zt 



where 



V{zt) = -Pexp 



■foo I 

00 



(8) 



(9) 



and ta are in the fundamental representation. Using (P) in the scattering amplitude 



gives 



{Q{Q)out\q{p)in) = {2n)6{p -q )u{q)-f u{p) / d^Zt 



V{zt) 



1 



J{qt~Pt)zt 



(10) 



The presence of the delta function in the amplitude is due to the target being (light- 
cone) time independent which leads to conservation of the "minus" component of 
momenta. It can be factored out in the standard fashion, 



{q{q)out\q{p)in) = (27r)5(p -q)M{p,q), 
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which gives the cross section 



da 



d^q 
(2^ 



(27r)5(2g+g- - q^,)e{q' 



2p' 



i27r)6{p--q-)\M{p,q)\' 



;i2) 



The local density of color charge in the nucleus, pa{xt,x~), is a stochastic variable 
which has to be averaged over. Commonly, one assumes a distribution of the charge 
sources which is local and Gaussian [Tl|, [16|, |1^], such that the average of any 
operator O is 



(0) = |l)pO[p]exp(-trpV/^^ 



(13) 



/i^(zt, X ) denotes the color charge density per unit transverse area d'^zt, and per unit 
of rapidity, dx~ /x~ , in the nucleus. When computing the total cross section, we will 



have to square the amplitude (10) before averaging over the color charge density p of 
the classical background field. On the other hand, for the case of elastic scattering, 
we first average the amplitude (p!0[ ) over p and then square it [1^, |I9[ . That way no 
color exchange occurs over a large distance in rapidity (from the projectile quark to 
the nucleus). 

The averages of V{zt) and V''{zt)V{zt) are given by |n 



{V{zt))p = exp 

and 

{V^{zt)V{zt)), = exp 



9\Nl - 1) , 

X / d yt G^[zt - yt) 



^'^^^ ^\ I d\ [Goiz, - y,) - Goizt - yt)? 



(14) 



■ (15) 



We have defined x{^' 



j^a/^o ^2(^^ -j^ which is the density of color charge 



in the nucleus per unit transverse area integrated over longitudinal phase space 
(rapidity), <^ Xq is the coordinate of the nucleus, with y^ = logXg/x^ its 
rapidity; while x^ <^ Xq is the coordinate of the quark projectile, and y = logx^/xg 
is its rapidity. ( reference point, which we choose to be midrapidity) . Go{zt — 

yt) is the free propagator of static gluons 



Go{zt - yt) 



A2 



d^kt e*^'(^'-s'*) 
(27r)2 kf 



(16) 
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The trace of the quark spinors in the squared amphtude is 



I 2 



(17) 



spins 



Averaging over the colors of the incoming quark is made trivial by the fact that (|T^JT^) 
are diagonal in color space. For elastic scattering, we shall use to color average 



the amplitude given by (10), and afterwards square it. The color averaging of the 
amphtude leads to a delta-function of transverse momenta (27r)^5^(pt — qt). This 
is due to the assumed translational invariance of the target in transverse space (we 
have assumed a large cylindrical target nucleus so that the color charge density is 
uniform) and should be understood as nR\{27i)^6'^{pt — qt) in the squared amplitude. 
Using ([T0|), ([n|) and ([ID finally leads to 



(Ph 



l-e 



QCD 



(18) 



where we have introduced the saturation scale of the target, = {N^ — l)alx/'^ 
(this is the same definition as in []TU|)- 

To compute the total cross section for scattering of the quark on the nucleus, 
we first square the amplitude (jT^) and then average over the background field us- 



ing (p^plSl) . It leads to 
d^qt 



d TtC 



2ttQ 



1^ Jd'^kt/kf [l-cxp(iktrt)] 



2^-n^QVN,AlcD + 1 



.(19) 



The integral over qt just gives S'^{rt) which in turn can be used to perform the 
integration. The final result for the total cross section is 



da^ 
d% 



1 _ g-'^'Q'/^-^QCD 



From eqs. ([T8|) and (2C), we see that in the high energy limit {Qs 



a 



tot 
qA 



2a\ 



el 
qA 



27iR\ 



(20) 

oo) we have 
(21) 



as a consequence of unitarity. 

It is more interesting to consider the differential cross section da^°^/ d%d'^qt when 
the quark momentum qt is large {qf ^ ^^qcd)- ^^is limit, one can convolute the 
quark-nucleus cross section with the quark distribution function in a proton and the 
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quark fragmentation function into hadrons, thereby relating qA scattering to single 
inclusive hadron (jet) production in pA collisions |^ . 



At very large transverse momentum, ^ Q^, we can expand the V^s in eq. (|15[), 
keeping only the first non-trivial term. This corresponds to the dilute perturbative 
limit. Using ( [16|) then gives the expected 

d^b ~ ^ • ^^^^ 

In the region q^ ^ ^qcd> turn, we must resum higher twists, i.e. rescatter- 

ings in the nuclear field. We then obtain 



dag A 1 

r\-/ 

d^qtd^h q^ 



(23) 



Clearly, the non-linearities of the classical field have flattened the differential cross 
section as compared to the perturbative cross section (22). This arises, in fact. 



from a suppression of particle production. It should be easy to distinguish from 
Prbroadening, i.e. "initial-state" interactions of the beam quarks with spectators 
from the target pO], which enhance the cross section at semi-high pf. 



The single inclusive p + A h + X cross section can now be obtained in principle 
by convoluting the total qA cross section with the quark distribution function of the 
proton at the factorization scale Q'j, which can for example be chosen to be qf, and 
the quark-hadron (jet) fragmentation function Dq/h{z,Qjr), where z is the ratio of 
hadron and quark momenta. 

J pA-*hX . J <7(T*°* 

Here, x denotes the fractional (light-cone) momentum carried by the quark, such 
that the light-cone coordinate of the incident proton is related to that of its quark 
by x~ = xx~ . 

The qA cross sections depend on x through the saturation momentum of the 
nucleus. From the analysis of HERA DIS data, the saturation momentum scales like 
Q'l{x~ / Xq) / QI{xq) = {xo/x~y = {x^x/x~y, with A ~ 0.3 For gold nuclei, 
it has been estimated that QHxq) ^ 2 GeV^ at RHIC energy (lOOA GeV gold 
beam, corresponding to ua = 5.4) [0. Therefore, near the rapidity of the incident 
proton, i.e. Up = loga;~/a;o = —5.4, we estimate that roughly ^ 10 GeV^. With 
1/Aq(^£, ~ 25 GeV~^, the exponent in eqs. (p!^,pO|) is < —800. In fact, even at 
6 ^ 6 fm from the center of a gold nucleus, using a Gaussian density distribution 
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for the nucleus we estimate that Ql{xQ,b ~ 6 fm) ~ 0.75 GeV^, and therefore 
Ql{x~,h ~ 6 fm) ^ 3.8 GeV^ ~ lOOAg^^,. At the edge of the acceptance of the 
BRAHMS spectrometer at RHIC, y = logx^ /x~ = 4, one obtains Q^{x~,b ~ 0) ~ 
6.6 GeV^. That is, 2(7^^ = cr*'^ = 2'7iR\, with Ra the radius of the black disc formed 
by the nucleus (i.e. the distance from the center of the nucleus where becomes 
of order Agf^j^), is quite close to the geometric cross section of the gold nucleus. 

The quark distribution functions in the proton, and their fragmentation into 
pions will modify the transverse momentum dependence of the pion cross section (p^ 
relative to the quark cross section (pSl) . Nevertheless, that modification is the same 
for both the dilute perturbative estimate of the qA cross section (^) as well as 
for the resummed qA cross section from (|23|) . Therefore, a modification of the 
Pt distribution from strong nonlinear color fields as compared to the dilute limit 
should hold independently of any scale dependence of the quark distribution and 
fragmentation functions. The big advantage of measurements in the forward region 
is that these effects extend much farther in pt than in the central rapidity region, 
hopefully making it much less ambiguous to observe them experimentally. 
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